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The CERN Axion Solar Telescope (CAST) searches for a —> 7 conversion in the 9 T magnetic 
field of a refurbished LHC test magnet that can be directed toward the Sun. Two parallel magnet 
bores can be filled with helium of adjustable pressure to match the X-ray refractive mass m 7 to the 
axion search mass m a . After the vacuum phase (2003-2004), which is optimal for m a < 0.02 eV, we 
used 4 Hc in 2005-2007 to cover the mass range of 0.02-0.39 eV and 3 He in 2009-2011 to scan from 
0.39-1.17 eV. After improving the detectors and shielding, we returned to 4 He in 2012 to investigate 
a narrow m a range around 0.2 eV (“candidate setting” of our earlier search) and 0.39-0.42 eV, the 
upper axion mass range reachable with 4 He, to “cross the axion line” for the KSVZ model. We 
have improved the limit on the axion-photon coupling to g a 7 < 1.47 x ICC 10 GeV -1 (95% C.L.), 
depending on the pressure settings. Since 2013, we have returned to vacuum and aim for a significant 
increase in sensitivity. 

PACS numbers: 95.35.+d, 14.80.Mz, 07.85.Nc, 84.71.Ba 


Introduction —The low-energy frontier of elementary 
particle physics includes numerous experimental ef¬ 
forts, ranging from the search for neutrino-less double 
beta decay and searches for nucleon and electron electric 
dipole moments all the way to searches for new low-mass 
bosons. The best-motivated case for the latter remains 
the axion, the pseudo Nambu-Goldstone boson of a new 
broken global 17(1) symmetry, that is required in the con¬ 


text of the Peccei-Quinn mechanism to explain why CP- 
violating effects are extremely small or absent in QCD 
[ 7 ), :8|. The axion is also an excellent candidate for the cold 
dark matter of the universe [§} and ongoing experimental 
searches as well as new efforts have recently gained fresh 
momentum U HI • On the other hand, the phenomenol¬ 
ogy of axions has inspired the proposition of axion-like 
particles (ALPs) and other WISPs (weakly interacting 



2 



m axio„( eV ) 


FIG. 1. Exclusion regions in the m a -g al -plane achieved by 
CAST in the vacuum phase [32;, 33] and with 4 He [34[ and 
3 He [35l . 36 ] filling. We also show constraints from horizon¬ 
tal branch (HB) stars [4ll . [42] (a similar limit stems from 
blue-loop suppression in massive stars [43|), and the hot dark 
matter (HDM) bound [US3- The yellow band represents 
typical theoretical models with \E/N — 1.95| = 0.07-7. The 
green solid line corresponds to E/N = 0 as in the KSVZ 
model [5(J, [5j|, a typical hadronic axion model. In red we 
show our new limits near m a = 0.2 and 0.4 eV from our 2012 
data taking campaign with 4 He gas. 


sub-eV particles) with additional theoretical motivations 
from string theory and cosmology. Experimental searches 
for low-mass bosons include precision searches for new 
long-range forces 2 # 27 and oscillating nucleon electric 
dipole moments j28H30( . 

The most promising technique to search for axions to 
date remains Sikivie’s idea to use the two-photon ver¬ 
tex of axions or ALPs and look for their conversion to 
photons in a strong external magnetic field B (also re¬ 
ferred to as the Primakoff effect) 1 3j]] - The smallness 
of the coupling strength g aj can be overcome by coher¬ 
ent conversion in a macroscopic volume. For more than 
ten years, the CERN Axion Solar Telescope (CAST) has 
pur sued this idea in the form of a large axion helioscope 
|32 - l36| . Earlier helioscope searches were conducted in 
Brookhaven 371 and Tokyo [38l-|40j. In CAST, a refur¬ 
bished LHC dipole magnet (length 9.26 m, field strength 
9.0 T) is mounted to follow the Sun for approximately 
1.5 h at both sunrise and sunset. In this way various 
X-ray detectors at both ends of the magnet are used to 
look at the Sun. 

For a fixed axion-photon coupling strength p a7 , the 
a -7 conversion probability decreases when the a -7 oscil¬ 


lation length becomes smaller than the magnet length, 
which limits CAST’s sensitivity to an axion mass range 
m a 0.02 eV. For larger m a values, the transition rate 
can be restored by providing photons with a refractive 
mass using a low-Z gas filling [52|. The gas density (i.e. 
pressure) defines the search mass and one needs to step 
through many pressure settings to search a broad m a 
range. 

CAST has taken data since 2003 and produced the ex¬ 
clusion plot shown in Fig. [T] The vacuum phase (2003- 
2004) achieved the best sensitivity up to m a ~ 0.02 eV 
because of the long exposure time. With 4 He filling 
(2005-2007), the maximum search mass was extended 
to ~ 0.4 eV, corresponding to the 4 He vapor pressure 
at cryogenic temperatures. With 3 He in 2009-2011 we 
finally reached 1.17 eV, which exceeds the cosmological 
hot dark matter limit for axions 44-49J. For m a >0.6 eV 
our 3 He search has reached the “axion line” for a typi¬ 
cal hadronic axion model [H ED] (E/N =0, see Fig. [TJ , 
i.e., the locus of g aj vs. m a which is motivated by QCD 
axions, as opposed to more general ALPs. 

To make the experimental progress pushing the sen¬ 
sitivity to lower coupling constants of Fig. [Q a much 
larger helioscope, the International Axion Observatory 
(IAXO) [531 - 1551] . has been envisioned and, for very small 
ro„,, a large-scale photon-regeneration experiment ALPS- 
II j56[ has been put forward. Until these next-generation 
projects become operational, CAST can push its bound¬ 
aries by reducing background levels and increasing the 
exposure time. Therefore, CAST has engaged in a pro¬ 
gramme of technical improvements and renewed data 
taking to achieve these goals. 

We report here on results from our 2012 data taking 
campaign with 4 He gas. At each pressure setting, we 
integrated for about 7.5 h per detector (5 solar track¬ 
ings), in contrast to 1.5 h (1 solar tracking) in our earlier 
4 He campaign. On the sunset side, new Micromegas de¬ 
tectors with improved shielding and veto were installed, 
decreasing the background level by an approximate fac¬ 
tor of 4 |57[ . The search masses were selected to cover a 
previous “candidate setting” ( m a ~ 0.2 eV), where un¬ 
usually many events above background had appeared, 
and the upper end of what can be reached with 4 He 
(m 0 ~ 0.4 eV) in order to reach the KSVZ (Kim, Shif- 
man, Vainshtein, Zakharov model) axion line. 

System description —After completing the 3 He phase, 
the gas was removed from the system. The magnet bores 
were then filled with 4 He to take advantage of the sophis¬ 
ticated gas metering system capable of filling the cold 
bore in small steps with a reproducibility of better than 
100 ppm [ 35 I ] . 

As explained in Ref. [36|, at increasingly higher gas 
densities (14-108 mbar at 1.8 K), 3 He gas dynamics (con¬ 
vection and buoyancy) at the ends of the cold bore affect 
the density distribution along the cold bore. This effect 
progressively shortens the region with uniform density 
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(effective coherence length) and also causes variations in 
the central density during the tilting of the magnet. In 
the 4 He case, the densities involved are relatively low (up 
to 15.5 mbar at 1.8 K) and these effects are negligible. 

For 4 He, the central gas density was calculated from 
the equation of state of the 4 He gas using the measured 
cold bore temperature, gas pressure and magnet vertical 
angle. The coherence length was taken as the magnet 
length of 9.26 m. 

Whilst Computational Fluid Dynamics (CFD) simu¬ 
lations have not been performed for 4 He, the 3 He CFD 
simulations at a similar pressure over various vertical an¬ 
gles gave effective coherence lengths above L e ff = 8 m. 
The main source of systematic error in the present result, 


according to our previous studies [33j, |34j, |36| is given by 
the uncertainties on the gas dynamics inside the magnet 
bores. In order to estimate the effect of such uncertain¬ 
ties the upper limit calculation has been redone using a 
reduced effective coherence length L e g = 8 m for all ver¬ 
tical angles. This value represents a worst case scenario 
of the effect of gas dynamics on L e g according to CFD 
simulations performed in Ref. ]36j. 

The X-ray detectors installed at CAST during the 2012 
data taking campai gn w ere three Micromegas detectors of 
the microbulk type 5^- 621 (one in the sunrise and two in 
the sunset side) and a pn-CCD detector in the focal plane 
of an X-ray telescope [63J on the sunrise side. While the 
detectors on the sunrise side remained unchanged since 
the previous data taking campaign, the Micromegas de¬ 
tectors on the sunset side were upgraded, improving the 
background levels of the detectors. This is the result of 
low background techniques developed for the Micromegas 
detectors [57|, where different strategies were exploited: 
the manufacturing technology of the novel microbulk Mi¬ 
cromegas, the intrinsic radiopurity of the detectors fbl- . 
the discrimination algorithms on the analysis, and differ¬ 
ent shielding strategies. 

The upgrade focuses on reducing the contribution of 
the environmental gamma flux. The lead shielding thick¬ 
ness was increased (from 25 mm to 100 mm) and the 
design is more compact, improving the shielding around 
the pipes to the magnet. The inner copper shielding was 
increased from 5 mm to 10 mm of copper. It is con¬ 
nected to the magnet bores by a 10 mm thick copper 
pipe, which has an inner polytetrafluoroethylene coating 
of 2.5 mm thickness in order to attenuate the 8 keV cop¬ 
per fluorescence. In addition, the aluminum strongback 
was replaced by a more radiopure copper one, and all 
the components close to the detector have been carefully 
selected and cleaned. A plastic scintillator was installed 
on the top of the shielding, allowing the discrimination 
of background events induced by cosmic muons in the 
detectors (see Ref. [5?]] for more details). 

These upgrades reduce the background level to 1.3 and 
1.7 x 10 -6 c cm' 2 keV -1 s _1 for the sunset detectors 1 
and 2, respectively, i.e., about a factor 4 lower than be- 



FIG. 2. Comparison between background (blue bars) and 
tracking (red bars) spectra of one Micromegas detectors in¬ 
stalled in the sunset side. 


fore, improving the signal-to-noise ratio by a factor of 2. 
Background and tracking spectra of one of the detectors 
in the sunset side are shown in Fig. [2] As seen both lev¬ 
els are compatible within their error. The background is 
dominated by the copper fluorescence at 8 keV, its escape 
peak above 5 keV and the Argon fluorescence around 3 
keV. According to our current understanding of the back¬ 
ground [57], [58|, [67| , these fluorescences are induced by 
secondary particles generated in the inner materials by 
external radiation, mainly muons that are not tagged by 
the active veto system. On the other hand, the detector 
response has been fully characterized by a detailed sim¬ 
ulation [67] and calibrations at different energies in an 
X-ray beam [66j. Moreover, daily calibrations allow us 
to monitor the performance of the detectors during the 
data taking period (see Ref. [57| for more details). 

The extensive research program for the background re¬ 
duction continues; different measurements in special set¬ 
ups have been performed, such as the ones underground 
at the Laboratorio Subterraneo de Canfranc 651, where 


the cosmic muons are highly suppressed, and the test at 
surface level These special test benches, together 
with the simulations [67j , were crucial in order to under¬ 
stand the origin of the background in the Micromegas 
detectors. Additional upgrades have been implemented 
at CAST for the re-scanned vacuum phase that started 
in 2013, like the AFTER 681 front end electronics and a 


new design of the active shielding on the sunset side, with 
improved background levels on the Micromegas detectors 
further. 

Data analysis and results .—The results presented in 
this work are related to the data acquired with three Mi¬ 
cromegas detectors during 2012 using 4 He as a buffer gas. 
The total exposure time in axion-sensitive conditions (so¬ 
lar tracking) was about 147 h per detector with a back¬ 
ground time of 2277 h per detector. In the first part 
of the 2012 data taking campaign, the axion mass range 
0.39 < m a < 0.42 eV was scanned, which corresponds 
to the pressure range 13.9-15.5 mbar at 1.8 K. During 
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this period, 17 pressure settings were covered with a step 
size of 0.1 mbar and an effective average exposure time of 
~24 h per setting. In the second part of the data taking 
campaign, a single setting at m a ~ 0.20 eV was covered 
during ~ 30 h. 

The analysis of the data was performed in the same 
manner as in our previous publications [35], [36J]. In order 
to take into account the variations of the density inside 
the magnet bores during the trackings, an unbinned like¬ 
lihood was implemented, where the likelihood function 
can be expressed as 


log£ oc -R t + ^2 log R(U, Ei , di ). 


( 1 ) 


Here, Rt is the expected number of counts from the axion 
to photon conversion over all the exposure time, energy 
and detectors. The sum is over each of the n detected 
counts during the tracking time, for an expected rate 
R(ti, Ei,di) as a function of the event time ij, energy Ei 
and detector di, given by the expression 


R(t,E,d) = B d + S(t,E,d), (2) 


where Bd is the background level of the detector d. 
S(t,E,d) is the expected rate from axion conversion in 
the detector d given by 


S(t,E,d) = ^Pa^e d . 


( 3 ) 


Here, e d is the detector efficiency, is the differential 
solar axion flux, which can be parametrized 41] by the 
expression 


dd>g 

dE 


6.02 x 10 lo gl o 


£2.481 

e S/1.205 


[cm 2 s 1 keV 


( 4 ) 


with g w = g a7 /(10 _lo GeV _1 ) and energies in keV. P 0 _> 7 
is the axion to photon conversion probability inside a 
strong magnetic field ]52|, given by 


Ra—t 7 — 


j„ 7 B\ 2 1 + e TL — 2e ri / 2 cos(gL) 
2 ) q 2 + T 2 /4 


( 5 ) 


where q = |m 2 — to 2 |/(2.E) is the axion-photon momen¬ 
tum transfer in a magnetic field B and T is the absorption 
coefficient in the buffer gas. 

As explained in Ref. [34] , the dependence on m a in the 
expressions above is included in the conversion probabil¬ 
ity given in Eq. ©, where P Q _>, 7 is enhanced for axion 
masses that match the refractive photon mass m 7 as de¬ 
termined by the buffer gas density. Therefore for a given 
axion mass m a , only the counts for which the coherence 
condition is fulfilled will contribute to the likelihood func¬ 
tion log/1. 

By maximizing log C a best-fit value gP n is obtained. 
This value is compatible with the absence of a signal in 




FIG. 3. Expanded view of the limit achieved during the 2012 
CAST data taking campaign, using 4 He as buffer gas. The 
plot on the bottom corresponds to the “candidate setting” at 
0.2 eV, while the plot on the top is the excluded region above 
0.4 eV. The green line represents the KSVZ benchmark model 
with E/N = 0. 


the entire axion mass range and thus an upper limit on 
g ai is extracted, by integrating the Bayesian posterior 
probability from zero up to 95% with a flat prior in g 4 7 - 
The computed upper limit for several values of m a is 
displayed in red in Fig. [I] A close view of the excluded 
region is shown in Fig. [3j where only the axion mass 
range scanned during the 2012 data taking campaign is 
included. 

This demonstrates that CAST improves its previous 
limit for axion masses 0.39 < m a < 0.42 eV, down to 
an average value of the coupling constant g ai ~ 1.47 x 
10 -10 GeV -1 . Moreover, the “candidate setting” around 
0.2 eV was rescanned and in the absence of an axion 
signature a coupling limit of g ai < 1.40 x 10 -10 GeV -1 
at a 95% of C.L. was derived. In contrast to the previous 
results the structure of the excluded region has 

a smooth shape, because of considerably larger exposure 
time per step, resulting in smaller statistical fluctuations. 

Conclusions- After finishing its initial mission of 
searching for solar axions up to the largest m a value that 
could be reasonably reached with 3 He, CAST has em¬ 
barked on a programme of technical improvements aimed 
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at increasing its sensitivity. Here we reported the first 
results from a new 4 He run with significantly reduced 
background rates in two narrow m a bands around 0.2 
and 0.4 eV. The bounds on g ai were significantly im¬ 
proved in these regions, allowing us to cross the KSVZ 
axion line near the highest possible 4 He pressure. 

Since 2013, we have returned to vacuum measure¬ 
ments, i.e., the low-mass regime, m a < 0.02 eV. With 
various further improvements, notably an additional X- 
ray telescope, we aim at a sensitivity for g ai — 0.6 x 
10” 10 GeV -1 for these low masses that could even sur¬ 
pass a recently improved stellar-evolution bound from 
the helium-burning lifetime of globular-cluster stars [42| ■ 
This low-mass regime, while far away from the “axion 
line,” is important in the context of ALPs, e.g. the prop¬ 
agation of TeV gamma-rays in large-scale astrophysical 
magnetic fields can be addressed by a possible photon- 
ALP conversion. 
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